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Abstract 

A systematic analysis of transverse momentum and rapidity distributions measured 
in high-energy proton - proton (pp) collisions for energies ranging from 53 GeV to 7 TeV 
using Tsallis thermodynamics is presented. The excellent description of all transverse 
momentum spectra obtained in earlier analyses is confirmed and extended. All energies 
can be described by a single Tsallis temperature of 68 ± 5 MeV at all beam energies 
and particle types investigated (43 in total). The value of the entropic index, q, shows a 
wider spread but is always close to q 1.146. These values are then used to describe the 
rapidity distributions using a superposition of two Tsallis fireballs along the rapidity axis. 

It is concluded that the hadronic system created in high-energy p - p collisions between 
53 GeV and 7 TeV can be seen as obeying Tsallis thermodynamics. 


1 Introduction 


Recently a power-law function based on the Tsallis distribution has been used extensively in 
hts of the pt distributions measured in high-energy collisions The relationship with an 

approach based on Tsallis thermodynamics |6]-[^ (and not simply the distribution associated 
with it) has been clarihed in several papers p|-[I^, in particular, the exponent appearing in 
the probability distribution used in 011 is associated with the entropic index, q forming the 


basis of Tsallis thermodynamics. A further step was made in 12 to generalize Hagedorn’s 
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theory 13,14 to non-extensive thermodynamics that predicts a limiting temperature (T) and 


a characteristic entropic index for hadronic systems 12 


In this paper a systematic analysis of transverse momentum and rapidity distributions mea¬ 
sured in high-energy proton - proton (pp) collisions for energies ranging from 53 GeV to 7 
TeV using Tsallis thermodynamics is presented. The excellent description of all transverse mo¬ 
mentum spectra obtained in earlier analyses is confirmed and extended. The present analysis 


of pT-distributions confirms previous results 10 ,11,15 ,16 giving constant temperature and 


constant entropic index when the non-extensive distribution for systematic analysis of experi¬ 


mental data. A recent analysis 17,18 obtained a higher temperature but it can be shown that 


this is due to not using a summation over all charged particles produced but simply fitting a 
single formula, doing so makes it consistent with the temperature T used in this paper. In this 
paper these values are used to analyse rapidity distributions. In this case we use a model to 
describe the fireball rapidity distribution along the rapidity axis and show that a consistent 
pattern emerges which allow us to make extrapolations for the upcoming pp collision at the 
upgraded LHC. 

In what follows we present a new analysis of py-distribution extending those already pub¬ 
lished, and provide more precise values for T and q. Then we present a model based on the 
theory for describing rapidity distribution and perform a systematic analysis of experimental 
data on rapidity distribution of charged hadrons is performed. Finally the results of the analysis 
are used for predicting and rapidity distributions for the next LHC phase at 13 TeV. 

2 Transverse Momentum Distributions 


The power-law distribution based on Tsallis thermodynamics is given, in terms of transverse 


momentum and rapidity, by 10 
(fN 


Pt dpT dy 


= gV 


rriT cosh y 
(2vr)2 


l + (g-l) 


rriT cosh y — p 
T 


q 

q-l 


( 1 ) 


where p is the chemical potential, rriT = \/pt + ^ volume, g is the degeneracy factor 

and q is the entropic factor, which measures the non-additivity of the entropy. Boltzmann 
statistics is recovered with q = 1. This equation can be obtained also from a non-extensive 


version of the perfect gas partition function 19 


It has been shown that at central rapidity y = 0 one can easily obtain the transverse 
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momentum distribution in terms of the central rapidity density, ^ 


, as 

20 

y=0 



(PN 


dprdy 


Pt rriT dN 


y=0 


T dy 


(2-g)(3-2g) 


v=o (2 - 0)^1 + 2moT + 2T2 _ 




9-1 


1 + (g- 1) — 


q 

9-1 


( 2 ) 


The above equation will be used for a systematic analysis of pr-distribution for pp collisions 
for a large set of experiments, as listed in Table A subset of those experiments was already 
analysed in Ref. 21 . Since all experiments report results for a narrow range of rapidity around 
the central region the approximation used is appropriate. The only free parameters to be 
adjusted to the experimental data are q and T. According to the thermodynamical theory in 


Ref. [^, if both the self-consistent principle from Hagedorn 13 and the Tsallis statistics [& 
can be applied in high-energy physics then both T and q must be independent of particle type. 


Table 1: 


Index 

Particle 

\/s(TeR) Exp. Ref. 

Index 

Particle 

^siTeV) 

Exp. 

Ref 

1 

7r+ 

0.2 

PHENIX 

22 


22 

P- 

0.9 

CMS 


24 

2 

7r“ 

0.2 

PHENIX 

22 


23 

A 

0.9 

CMS 


5 ; 

3 

K+ 

0.2 

PHENIX 

22 


24 


0.9 

CMS 


5 

4 

K- 

0.2 

PHENIX 

22 


25 

'7r+ 

2.76 

CMS 


24 

5 

P+ 

0.2 

PHENIX 

22 


26 

'K~ 

2.76 

CMS 


24 

6 

P- 

0.2 

PHENIX 

22 


27 

K+ 

2.76 

CMS 


24 

7 

P+ 

0.2 

PHENIX 

22 


28 

K- 

2.76 

CMS 


24 

8 

P- 

0.2 

PHENIX 

22 


29 

P+ 

2.76 

CMS 


24 

9 

7r° 

0.9 

ALICE 

23 


30 

P- 

2.76 

CMS 


24 

10 

'7r+ 

0.9 

ALICE 

3 


31 

7r° 

7.0 

ALICE 


23 

11 

'K~ 

0.9 

ALICE 



32 

'7r+ 

7.0 

CMS 


24 

12 

vr^ 

0.9 

CMS 

24 


33 

TX~ 

7.0 

CMS 


24 

13 

TX~ 

0.9 

CMS 

24 


34 

K^s 

7.0 

CMS 


5 ; 

14 

K 

0.9 

CMS 

5 


35 

K+ 

7.0 

CMS 


24 

15 

K+ 

0.9 

ALICE 

3 


36 

R- 

7.0 

CMS 


24 

16 

K- 

0.9 

ALICE 

3 


37 

V 

7.0 

ALICE 


23 

17 

K+ 

0.9 

CMS 

24 


38 

K* 

7.0 

ALICE 


25 

18 

K- 

0.9 

CMS 

24 


39 

P+ 

7.0 

CMS 


24 

19 

P+ 

0.9 

ALICE 

3 


40 

P- 

7.0 

CMS 


24 

20 

P- 

0.9 

ALICE 



41 

0 

7.0 

ALICE 


25 

21 

P+ 

0.9 

CMS 

24 


42 

A 

7.0 

CMS 
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Figure 1: Fits to the pT-distribution given by Eq. [^to experimental data, (a) 7r+, ^/s = 0.2 


TeV, (b) A, ^/s = 0.9 TeV, (c) iF , ^/s = 2.76 TeV, 24 and (d) 0, ^/s = 7.0 TeV, 25 


In Fig. [^we show a typical result for the fitting of Eq. [^to measured pT-distribution. With 
the fittings for all experimental data used here we obtain the values for T and q that are plotted 
in Fig. 1^ where we observe a distribution of both temperature and entropic index around a 
constant value. 

The results obtained here are in agreement with previous analysis both qualitatively and 
quantitatively. Therefore we can determine the mean values from the analyses [lO|[T5|[l6|[^ , 
obtaining q = 1.146 ± 0.004 and T = 68 ± 5 MeV. With these values all the thermodynamical 
aspects of the hot fireball are determined. 
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Figure 2: Results for the temperature, T and for the entropic index, g, for all the sets of data 
analysed in this work. The numbering on the abscissa corresponds to the entries in Table 1. 
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3 Rapidity Distribution 


In the Bjorken model 26 a plateau in the central region of the rapidity distribution exists 


reflecting the invariance of the fireball under Lorentz transformation for boosts in the beam 
direction. Experimental results for high energy collisions, however, have shown that such a 
plateau does not exist, data showing two small peaks symmetrically located near the central 
rapidity region. 


A systematic analysis performed in Ref. 27 shows that color-glass-condensate (CGC) ap¬ 


proaches give a good description of the results at the central region but fail to reproduce the 
data at large values of \y\. Here we use a model for describing the rapidity distribution in pp col¬ 


lision using the non-extensive distribution. A similar method was already used in Refs 28,29 


in terms of Boltzmann-Gibbs statistics but the present analysis keeps some differences with 
respect to the previous one, as will be detailed below. 

Suppose that the fireball is moving in the laboratory frame in the beam direction with 
rapidity i/f. According to Eq. [^the yield of secondaries from the decay of this hreball is given 
by 


d^N 


Pt dpT dy 


= gV 


rriT cosh y' 
(2vr)2 


l + (g-l) 


rriT cosh y' — p 
T 


q 

9-1 


( 3 ) 


where y = y' + yj. 

From now on we assume /i = 0 since we restrict ourselves to high energy pp collisions, and 
adopt a model for the hreball longitudinal expansion which is similar, but not identical, to 
the Bjorken scenario 26 . In the present model the system the basic hreballs have rapidity 


distributed according to the function z/(|//) so that 


d'^N 


N 

A 


= ^ / zz 


(Pf) X 


niT cosh(|/ -yf) 


1 + (g- 1) 


rriT cosh(?/ — yf) 


T 


9 

9-1 


dyf (4) 


Pt dpT dy A " (27r)2 

where i'{yf)dyf is the number of hreballs with rapidity between yj and yf + dyf , N is the fully 
integrated particle multiplicity and A is a normalization constant ensuring that 

1 /•“ d^N 

N 


-dpT dy = 1. 


1-00 Jo dPT dy 

It is possible to integrate Eq. |^over the transverse momentum pt, resulting in: 


( 5 ) 


1 dN 


1 

A 


'(Vf) 


sedi^iy -yf) 


2g-3 

{q- l)cosh(y-y/)\ 'j-i 


{mo{q -1) + Tsech{y - yf)) 


mo + 


T sech(y-yf) 

9-1 


q-l 


{-ml{q - 2) -h 2moT sech{y - yf) + 2T'^ sech^{y - yf)) 


47r2(g - 2){q - l)^{2q - 3) 


dyf. 


( 6 ) 
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The equation above gives the rapidity distribution of particles produced in pp collisions and 
is the basis for describing the available experimental results at different beam energies. 

For completeness we still need the distribution function v{yf)- Based on the shape of the 
experimental distributions we choose the following ansatz: 




( 7 ) 


where 


Gq{yo,a;yf) = 




TTCr 


(yf - yo? 

2^2 


( 8 ) 


and eq{x) is the q-exponential function dehned as 

eq(x) = [1 — (g — . 


(9) 

In Eq. [^i/o and a are respectively the peak position and the width of the q-Gaussian function and 
are considered here as free parameters to be adjusted in a systematic analysis of experimental 
data. Eq. [^assumes that the main hnal state in the collision is composed of two hreballs moving 
in the beam direction with opposite rapidities ?//, each of them being composed by hreballs 
with rapidity-gaps with respect to the corresponding cluster that are distributed according to 
a q-Gaussian function. 

Notice that q' in Eq. does not need to be equal to q and in the present work we perform 
the analysis of experimental data with two different assumptions: q' = q and g' = 1, the 
last assumption corresponding to a Gaussian distribution of the rapidity-gap of the hreballs in 
each of the moving cluster. In addition, we assume that the hreballs are described by the non- 
extensive self-consistent thermodynamics with the limiting temperature T and the characteristic 
entropic index, g, that were found in the analysis of pr-distributions and hadronic spectrum, as 
described above. The values adopted here are T =68 MeV and g =1.146. For simplicity we will 
adopt mo = 139.59 MeV, since most of the charged hadrons produced at high energy collisions 
are pions. The set of experimental data used for the rapidity analysis is shown in Table 

The main diherence between the analysis performed in the present work and that performed 
in Refs. 28,29 is, of course, the use of Tsallis statistics in the thermodynamical description 


of the hreball however other differences must be noticed. While in both temperature and 
chemical potential could depend on the hreball rapidity, here we adopt constant temperature at 
the limiting value predicted theoretically and null chemical potential regardless of the hreball 
rapidity. This approach is justihed by the fact that our analysis focuses on charged hadrons 
distributions produced at very high beam energies (> 500 GeV), which are dominated by pion 
production. This assumption can be improved with the availability of rapidity distributions for 


identihed hadrons 28 - 30 
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Table 2: 


i/i (GeV) Exp. Range Ref. 


53 

UA5 

y 

< 

3.1 

31 

200 

UA5 

y 

< 

4.6 

31 

200 

PHOBOS 

y 

< 

5.3 

32 

410 

PHOBOS 

y 

< 

5.3 

32 

546 

UA5 

y 

< 

4.8 

31 

630 

UA5 

y 

< 

5.5 

33 

900 

UA5 

y 

< 

4.6 

31 

900 

CMS 

y 

< 

2.5 

34 

1800 

CDF 

y 

< 

3.5 

35 

2360 

CMS 

y 

< 

2.5 

34 

7000 

CMS 

y 

< 

2.5 

36 



y 


Figure 3: Fit using Eq. [^to experimental data for rapidity distribution from pp collision at 900 
GeV from the UA5 collaboration 31 . Data for y < 0 are obtained symmetrically from |/ > 0. 
The full line represents the best fit using the q-Gaussian function (g' = q) and the dashed line 
corresponds to the best £t using the Gaussian function {q' = 1) in Eq. The dotted line shows 
the results when a delta function is used for the hreball rapidity distribution, i.e., viyf) = S{yf). 


In Fig. I^we show the results of hts to the rapidity distributions of particles produced in pp 
collisions from several experiments with collision energies ranging from 50 GeV up to 7 TeV. 
The set of experimental data used in the complete analysis is shown in Table We observe a 
good £t of formula to the experimental data in the whole range of rapidity available, which 
covers almost the entire rapidity region. These results and the others presented in Appendix A 
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shows that a consistent description of the experimental data can be obtained with the model 
described above. Notice that the only free parameters in the htting procednre are a and yo 
appearing in Eq. From the analysis performed on all the sets of experimental data in Table 
[^we get the valnes for those parameters that are presented in Fig. 

For both qnantities we observe a sharp increase np to collision energies aronnd 500 GeV. 
The nncertainties above this energy are also larger dne to the fact that the range of rapidity 
available is smaller. The striking featnre of these resnlts is the fact that both the peak position 
yo and the width a of the q-Ganssian remain within the same range i.e. are approximately 
constant for collision energies from 500 GeV to 7 TeV, while the beam rapidity varies by a 
factor 2 in this range of collision energies. These resnlts allow ns to conclnde that after the two 
original protons collide two clnsters of hreballs are formed and they move parallel to the beam 
direction with opposite rapidities given by |?/o| = 2.3. Each clnster is formed by hreballs that 
move with respect to the center of each clnster in the beam direction with rapidity distribntion 
that is described by a q-Ganssian fnnction with width a ~ 1.8. We will assnme in the rest of 
this paper that the valnes of yo and a will not change drastically for collision energies np the 
14 TeV. Of conrse, it is expected that the position of yo will increase at mnch higher beam 
energies. 

Below ~ 550 GeV the scaling observed above disappears, as can be seen in Fig. It is 
possible that at these collision energies the temperatnre achieved by the system is signihcantly 
lower than the limiting temperatnre, T, hence more stndies on the direction of verifying the 
role played e.g. by the chemical potential at low energies must be carried out, as performed 
in 


28 30 


When we use the Gaussian function instead of the q-Gaussian for the rapidity distribution 
of the hreballs, i.e., 

^(l//) = G{yo, Vf) + G{-yo, ct; yf) , (10) 


with G(|/o, cr-,yf) being the Gaussian function with peak position at yo and width a, the results 
of the htting procedure Eq. to the experimental data are very similar to those obtained 
with the q-Gaussian function, as shown in Fig. |^with dashed line. The best ht values for the 
parameters yo and a are also practically the same and they are presented in Fig. |^with open 
symbols. 


4 Energy Dependence of Multiplicity 

Once the parameters T, q, yo and a are determined, if we know the function N{^/s) giving 
the multiplicity as a function of energy we can determine the pt and rapidity-distributions 
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Figure 4: 


Values for the parameters a and yo in Eq. as a function of the collision energy. 


of secondaries produced in pp collisions at any collision energy sufficiently high to allow the 
assumption that those parameters are independent of beam energy and particle properties. 
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Therefore now we investigate how the function N{^/s) behaves. 

The energy carried out by all particles with mass m is determined by 

gV /■“ f°° (PN 

= / ^{PT,y)^ — —dpTdy, ( 11 ) 

[2TTy J-ooJo dprdy 

where 

^{Pt, y) = rriT cosh(j/) (12) 


The total energy carried by baryons is 


E = 


p{m)E{m)dm 


( 13 ) 


where p{m) is the baryon mass spectrum. The function p{m) was determined by the non- 


extensive self-consistent theory 12 and is given by 

{q — l)m 


p{m) = 


7 


m 


5/2 


1 + 


T 


1 

q-1 


(14) 


In Ref. 21 it was shown that this formula reproduces correctly the observed hadronic mass 
spectrum up to a mass of approximately 2 GeV. 

Substituting Eq. |^in Eqs. [^and IR we obtain a relation between the total baryonic energy 
and the multiplicity that is given by 


£ = /V(£)l r p{m) 


poo POO 


’ —oo J 0 


d?N 

Ppt, y) -J — 1-dpT dy dm . 
dprdy 


(15) 


Notice that for sufficiently high energy, when T, q, a and yo do not depend on the collision 
energy, the only term which can depend on the total energy is the multiplicity, N{E). At the 
limiting temperature the calculations are difficult to be performed because of the singularity 
in the partition function corresponding to the phase transition. Therefore we use a sum over a 
subset of the particles produced in high energy collisions, so that 


N{E) = KE, 


( 16 ) 


where 


K = 


A 


E 


POO POO 


i=i "'o 


d'^N 1 

£{PT,y) -1 - y- (mi)dpT dy 

dpT dy 


( 17 ) 


where the sum is over all the particles indicated in Table m* is the particle mass, and 


d^N , , 

is the double differential yield given by Eq. when the mass is used. 


( 18 ) 
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Table 3: 


Particle 

Mass (GeV/c^) 

TT^ 

0.140 


0.494 

P 

0.770 

K* 

0.892 

P 

0.938 

E+ 

1.189 

E- 

1.197 

A 

1.232 


1.321 

S* 

1.385 

TT* 

1.533 

n- 

1.672 


1.869 


1.971 

D* 

2.010 

K 

2.281 


5.271 


With Eq. IT and using all particles listed in Table I we get K = 0.009 GeV h This value 
is above the best htting value K = 0.004 GeV“^, showing that the more massive particles not 
included in our calculation give a signihcant contribution to the total energy. In Fig. we show 
the linear variation of multiplicity using the calculated value for K as & dashed line. 

In Fig. I^we present the multiplicity N obtained from the httings of Eq. [^to the rapidity 
distribution data using the q-Gaussian function for the hreball rapidity distribution (black 
symbols). We observe that after a fast increase for collision energies up to ~ 500 GeV the 
multiplicity increases more slower at higher energies. As already discussed, for low energies we 
cannot ensure that T is constant and that the chemical potential is null, therefore we will not 
deepen our discussion for the results at those energies. Nonetheless we htted the multiplicity 
data with two different htting formulas: one linear for energies above 500 GeV and the other 
logarithm for the full rapidity range. The results are plotted in Fig. 

We observe that the linear £t has a coefhcient very different from K calculated above. 
Indeed, the curve one would obtain with this coefhcient is the dotted line also shown in that 
hgure, which gives a consistent description of the multiplicity in the energy range from 500 GeV 
up to 2500 GeV but would not describe correctly the multiplicity at 7 TeV. This happens 
because the coefhcient K consider only the contribution from the hadrons listed in Table 
which results in an overestimation of the total energy carried by all hadrons according to Eq. [T^ 
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Figure 5: The particle multiplicity N as & function of ^/s. The solid line is a linear ht to the 
data for ^/s > 800 GeV. The dashed line represents the logarithm function. The dotted line is 
the multiplicity as a function of energy when the coefficient K is used. 

The problem is that as the collision energy increases, small variations in the temperature, 
which is tending to the limiting temperature, will result in a sharp increase of the number of 
massive hadrons, as described by p{m), and these particles contribute signihcantly to the total 
energy even if they have small contributions to the total multiplicity. Therefore the linear £t 
in Fig. I^can be estimated as an upper limit to the multiplicity for energies above 7 TeV. The 
logarithm £t, on the other hand, could be considered as the lower limit, since it systematically 
underestimate the charged hadrons multiplicities for collision energies above 2 TeV. 

As examples of application of the present model we now calculate the pt- and ^/-distributions 
for charged particles at collision energy E = 13 TeV, which is the expected energy at the LHC 
experiments. 

The rapidity distribution is shown in Fig. [fusing the linear extrapolation (full line), which 
gives a total multiplicity M = 83, and the logarithm extrapolation, which gives M = 55. These 
results are obtained with q-Gaussian functions for viyf) with //o = 2.4 and a = 1.8. 

The estimates for pr-distribution of charged particles are shown in Fig.[^ The lines indicate 
the lower and the upper limits of our theoretical calculations using Eq. and considering the 
uncertainties on the parameters q and T. The uncertainty in the multiplicity is not taken into 
account but we present the results obtained with both extrapolations for N{E), namely, the 
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Figure 6: The charged hadrons rapidity-distribution predicted for 13 GeV pp collisions. The 
full line corresponds to the multiplicity calculated through the linear extrapolation and the 
dashed line corresponds to the multiplicity obtained with the logarithm extrapolation, as given 

in Fig. 1^ 
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Figure 7: The charged hadrons p^-distribution predicted for 13 GeV pp collisions. The shaded 
regions show the expected distribution considering the uncertainties in the parameters T and 
q. The upper region with full lines corresponds to the multiplicity calculated through the linear 
extrapolation and the bottom region with dashed lines corresponds to the multiplicity obtained 
with the logarithm extrapolation, as given in Fig. 
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linear (full lines) and the logarithm (dashed lines). The shaded area between the lines show 
the region where we expect the data to be measured at LHC. 

5 Conclusions 

In this work we presented an analysis of the pt and rapidity distribution in pp collisions 
assuming that the hnal state can be described by two hreballs obeying Tsallis thermodynamics 
having a Tsallis temperature T = 68 ± 5 MeV and an entropic index q = 1.146 ± 0.004 
independent of beam energy. 

Our results on the transverse momentum distributions are in agreement with previous anal¬ 
ysis and allows a complete characterisation of the thermodynamics of hreballs. 

We developed a model to describe the rapidity distribution with non-extensive thermody¬ 
namics which is based on the assumption that the hot system formed after pp collisions can be 
represented as two clusters of hreballs moving in the beam direction with rapidities distributed 
according to a function vijjf). We use a sum of two q-Gaussian functions to describe 
It results that for energies above ~ 500 GeV all rapidity distributions can be reproduced with 
constant peak-position and constant width for the q-Gaussian functions. 

With these results we estimated the behaviour of the multiplicity as a function of the 
collision energy, allowing us to present estimates for the future LHG experiments at 13 TeV. 
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Figure 8: Fits to rapidity distributions measured in high-energy p-p collisions (see Table 2) with 
Eq. The experimental data are from: (a) UA5 collaboration (b) UA5 collaboration 


(c) PHOBOS collaboration 1^ (d) PHOBOS collaboration 1^ (e) UA5 collaboration 31 
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(e) 


Figure 9: Series of fittings of experimental data on rapidity distribution (see Table 2) with Eq. 

The experimental data are from: (a) UA5 collaboration (b) CMS collaboration (c) 
CDF collaboration (d) CMS collaboration (e) CMS collaboration 36 
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